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Abstract 

We reconsider the Hung and Shuryak arguments in favour of searching for the 
deconfinement phase transition in heavy ion cohisions downward from the nominal 
SPS energy, at Eiab ~ 30 GeV/A where the fireball lifetime is the longest one. 
Using the recent lattice QCD data and the mixed phase model, we show that the 
deconfinement transition might occur at the bombarding energies as low as Eiab = 
3 — 5 GeV/A. Attention is drawn to the study of the mixed phase of nuclear matter 
in the collision energy range Eiab = 2 — 10 GeV/A. 
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Over the last ten years, a fundamental QCD prediction of the phase transition from 
hadrons into a state of free quarks and gluons (quark-gluon plasma) has been studied ac- 
tively. Extensive lattice QCD calculations not only allowed to specify the deconfinement 
temperature, the order of phase transition and its flavor dependence both for pure gluon 
matter and for the plasma with dynamical quarks, but also gave insight into physical 
phenomena above and below the deconfinement temperature. By now, there is a rather 
long list of various signatures which can signal the quark-gluon plasma formation. Many 
of these signatures were tested experimentally at CERN-SPS with the ^^O and ^^S beams 
at 200 GeV/A where crucial conditions for the deconfinement transition are expected to 
be reached. Indeed, some predicted effects such as the strangeness enhancement, J/\l/ 
suppression or (f)/{p + uj) enhancement have been observed but their hadronic interpreta- 
tion cannot be excluded. General belief is that for more definite conclusions heavier ions 
and higher beam energies should be used. 

Recently, a new strategy of experimental search for the QCD phase transition in heavy 
ion collisions has been advocated The equation of state (EOS) is very 'soft' in a narrow 
range of temperatures (AT ~ 10 MeV) around the transition temperature, leading to a 
significant reduction in transverse expansion of the fireball formed in heavy ion collisions 
1^. This 'softness' of the EOS, estimated by Hung and Shuryak at the energy density 
Esp — 1.5 GeV/ fm^ and the beam energy Eiab ~ 30 GeV/A |jl|, affects not only the 
transverse but also the longitudinal expansion and results in a longest lifetime of the 
excited system. Radical changes in the hydrodynamic space-time evolution around the 
'softest point' of the EOS should lead for certain observables to a sharp and specific 
dependence on the heavy-ion beam energy . Therefore, some aspect of the deconfinement 
transition might be better studied at lower collision energy, downward from the nominal 
SPS energy. 

Hung and Shuryak are confident that even in more sophisticated models than used in 
ref.[|l|, the total lifetime of a fireball should have a maximum near the indicated collision 
energy region. From our point of view, this conclusion is not evident because of the 
two crucial assumptions made. Firstly, the crossover behaviour of the deconfinement 
transition was simulated in ref.|]I[ by an arbitrary smoothing of the results of a simple 
two-phase model exhibiting the first order phase transition. It is not clear how well this 
description approximates the lattice QCD results. Secondly, the EOS of the baryon-free 
matter, hb = 0, has been used |]l[. Due to a considerable stopping power, there is little 
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hope to create a baryonless matter in nucleus- nucleus collisions up to Eiab — 200 GeV/A 
(especially for heavy systems). In this paper, we shall consider how the longest-lived 
fireball is sensitive to the details of the EOS and how it survives in baryon-rich matter 
inherent in systems formed at energies Eiab < 200 GeV/A. Our analysis is based on the 
statistical model taking properly into account the mixed phase in which unbound quarks 
and gluons coexist with hadrons P, 

The EOS is the key quantity and thermodynamic properties of excited matter near 
QCD phase transition should be calculated from the first principles in a non-perturbative 
manner. At present, such calculations are possible only in the lattice QCD. This approach 
shows that gluonic matter in pure gauge SU{2) and SU{3) theories exhibits a phase 
transition of the second and first order, respectively (for a recent review see ref.|Q). For 
the more realistic case of the SU{3) theory with dynamical quarks, recent lattice studies 
IP show that there is a smooth crossover rather than a distinct phase transition if the 
quark masses are close to the physical ones. For the most important case of the lattice 
QCD at finite baryon number density, nothing is really known and statistical models 
should be invoked to describe thermodynamic properties of the excited nuclear matter. 
There are many versions of the statistical model (see ref . and references quoted therein) 
but all of them predict the deconfinement phase transition of the first order and, therefore, 
do not reproduce the lattice QCD results with dynamical quarks mentioned above. The 
only exception is the statistical mixed phase model of the deconfinement |^, which 
we will use in our considerations. The application of the mixed phase models to the 
deconfinement transition is supported strongly by the discussion of deconfinement as a 
color screening effect and the observation that the hadron-like excitations survive above 
the critical temperature (see for example ref. j^). 

The specific feature of the approach developed by Shanenko, Yukalova and Yukalov 
(SYY) 1^, ^ is to consider the coexistence of spatially non-separated the hadron phase 
and the quark-gluon plasma phase. The latter one consists of unbound 'generic' particles 
(quarks and gluons in our case) while hadron phase is composed of quark-gluon clusters. 
In the mixed phase model ||^, one assumes beforehand the separation of cluster degrees of 
freedom, i.e. an exact Hamiltonian H{ip) is replaced by an effective cluster Hamiltonian 
Hc{ipc)- Here, ip denotes field operators of the generic particles and ipc = {4'n', n = 1,2, . . .} 
stands for quasiparticle operators {n > 1 corresponds to clusters and n = 1 to unbound 
generic particles). A large variety of quasiparticles leads to an enormous number of 
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possible states. In the equilibrium, the physical state of a system corresponds to the 
extremum of the thermodynamic potential F{Hc). 

The following two important points should be emphasized. Firstly, the effective Hamil- 
tonian Hc{ipc) may acquire extra dependence on thermodynamic parameters like temper- 
ature T and cluster densities pc- He = Hc{ipc,T, pc) where pc = {pn]n = 1,2, . . .}. The 
appearance of a density-dependent interaction is a distinctive feature of the SYY approach 
as compared to other statistical models P Q Secondly, one must ensure that H{i/j) and 
Hciipc) are thermodynamically equivalent, i.e. their thermodynamic characteristics are 
the same in the thermodynamic limit of ^ oo keeping constant pc [§] . These demands 
of the thermodynamic equivalence and the thermodynamic self-consistency impose addi- 
tional conditions: 

hm i[F(if)-F(iJ,)]=0, 

V — >QO W 

PQ^const 

hm ^[dFiH)-dFiH,)] = (1) 

Pc — const 



which lead to 



dT 



,dpn, 



0, (2) 



n = l,2,... (3) 



and essentially define the form of the cluster Hamiltonian. In the mean-field approxima- 
tion, the cluster Hamiltonian becomes [HI: 



^c = EE /^r ^n(r,s) {Kn + Un{T,p,)) M^,s)-C-V , (4) 

n s 

where n enumerates the clusters and s stands for their internal degrees of freedom. Kn = 



— + is the kinetic energy and f/„ is a mean field acting on the n-particle cluster. 
In the same approximation, the condition of thermodynamic equivalence becomes: 

dUn{T,p, 



dT 



n = l,2,..., (5) 



^Similar situation is met in the case when interaction is taken into account by the excluded- volume 
method. 
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i.e. the mean fields Un = Un{pc) depend on the temperature only through the temperature 
dependence of the densities pn- The c-number term C - V ui is necessary to satisfy 

the conditions (0) and in which case C is a function of f/„'s [||. Unbound particles 
are treated as trivial clusters with n = 1. 

If masses M„ of the isolated clusters and their quantum numbers are known either 
experimentally or from other calculations, then, to apply the mixed phase model, one 
has to define only Un- The mean field Ui acting on unbound gluons or quarks can be 
approximated as follows 0, ^ : 

A 

Ui = — < 7 < 1 (6) 

where p = J2n ^ Pn is the total quark-gluon density. The presence of p in @ corresponds 
to the inclusion of the interaction between all components of the mixed phase. Note 
that when there are no hadrons, i.e. p„ = for n > 1, the expression (P) is just the 
same as used earlier to describe the thermodynamic properties of the quark-gluon plasma 
with a density-dependent quark mass For the case of n > 1, alongside with the 

mean-field term depending on the constants A and 7, there is also a cluster-cluster 
interaction potential which is proportional to the hadron density and the function $„m(^) 
characterizing the interaction strength between clusters of n and m generic particles ^ . 
For long-ranged cluster interactions it is possible to get a recurrent relation P| : 

$nm(r) ~ nm $o(r) , (7) 

and to reduce all unknown interactions to the single interaction potential ^o{r) between 
simplest non-trivial clusters, e.g. two-gluon glueballs in the ground state, lightest mesons 
or baryons. In the Hartree approximation to He, one needs to know only one constant: 
$0 = J dr $o(^) to describe the cluster-cluster interactions. Thus, the Hamiltonian (^) 
is completely defined and thereby any thermodynamic characteristics of the mixed phase 
system can be found if three parameters A, 7 and $0 have been fixed [0, This 3- 
parameter set was found 0] by fitting a temperature dependence of the energy density e 
and the pressure p calculated within the lattice QCD for the pure gluonic matter in the 
gauge SU{2) and S'f/(3)|T0[ theories. It is worth emphasizing that A depends only 
on the colour group and 7 is constant for all gauge systems. So, fitting the SU{3) pure 
gluonic QCD lattice data allows to fix A and 7 parameters which can then be used for 
the SU (3) system with quarks. As to $0 in this case, it can be found from a nucleon- 
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nucleoli potential by referring to the relation (|^). In the following, we use 7 = 0.62, 
^1/(37+1) = 225 MeV and $0 = 4.1 ■ 10"^ MeV-'^. 

The mixed phase model predictions for the SU (3) system with two light flavours and 
hb = ^ are shown in Fig. 1. When confronted with the lattice QCD data, the mixed 
phase model gives a very similar temperature dependence of the energy density e 
and pressure p. As follows from these results, the mixed phase model estimates the 
deconfinement temperature to be T^ec = 150 MeV and predicts the crossover-type phase 
transition P, Q in full agreement with the QCD lattice data. Some overshooting of the 
lattice data at T > 1.5 Tdec is related to neglecting the negative Coulomb-like term of 
the quark-gluon interactions. Two different sets of the best available lattice QCD results 
plotted in Fig. 1, correspond to two different scheme of including quarks. A peak-like 
structure of e near the deconfinement temperature T^ec for the Kogut-Susskind scheme 



13] seems to be unrealistic, vanishing for temperatures below ~ 0.9 T^ec- The quark mass 



in the Kogut-Susskind calculations amounts to 0.1 Tdec which is close to the physical 
mass used in the mixed phase model. In the Wilson scheme we do not really know the 
quark mass used. In this case, the value ~ T^ec given in [0 seems to be enormously 
large because for T < T^ec there is a good agreement between the Wilson scheme results 
and the ideal meson gas calculations |12| . 

The EOS in the form advocated by Hung and Shuryak |Q is represented in Fig. 2. 
Here, all data of Fig.l are replotted alongside with the curve used in to simulate 
the crossover transition. The p/e-f unctions for lattice data themselves show a minimum 
which is just associated with the softest point of the EOS where a fireball of the excited 
nuclear matter lives longest. In the p/e-representation, the lattice QCD data for two 
schemes of accounting for dynamical quarks are strikingly different but, nevertheless, 
they yield the same position of the minimum at Esp ~ 0.5 — 0.6 GeV/ fim? which can 
be reached with the Au beam energies of about 3 — 5 GeV/A. It is noteworthy that 
this value od Esp corresponds roughly to the value of the bag constant, so the fireball 
near the softest point may be considered as a 'big hadron'. The mixed phase model 
predicts the position of a minimum near Esp ~ 0.3 GeV/ fm? which is rather close to 



sp ~ 



the lattice results. On the contrary. Hung and Shuryak obtain much higher value, e 
1.5 GeV/ fm? |jl| and the minimum is seen more distinctly than in the lattice data. Since 
E ~ T^, the difference in deconfinement temperatures is small. The Kogut-Susskind and 



Wilson schemes predict the deconfinement temperatures T^ec as high as 157 MeV [13 
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and 150 MeV ^ respectively, and the mixed phase model gives Tdec = 150 MeV. The 
Hung-Shuryak approximation does not come from a fit to the lattice QCD data but is an 
arbitrary smoothing of the two-phase model results where T^ec = 160 MeV is defined by 
the choosen bag constant 0. 

The mixed phase model can be naturally generalized to the case of the non-zero baryon 
density ub [Q- The n^-dependent pressure in the SYY approach coincides with that of 
the Walecka-like model up to ra^ ^ 3 [n|. At higher baryon densities, the mixed phase 
model gives lower values of pressure due to the quark admixture. As to the predicted 
energy density, the difference between these two models is even smaller. 

As is seen in Fig.3, the p/e-function changes drastically with increasing the baryon 
density of a system: the position of Egp remains unchanged at ra^ < no though the 
minimum is gradually vanishing with increasing ub and disappears at ub ^ . So, 
we arrive at somewhat controversial demands: to reach the condition of the longest-lived 
fireball in heavy ion collisions one should go downward in bombarding energies as far as 
Eiab ~ 3 — 5 GeV/ A but, on the other hand, high baryon density of a fireball formed at 
these energies will suppress much the effect of the softest point. Due to smallness of Sgp 
it may be more favourable to look for the softest point effect in non-central high-energy 
collisions in the target fragmentation region. To see how the lifetime of a fireball will 
change with the beam energy, one needs detailed dynamical calculations with EOS of the 
mixed phase model Q. It is of interest to note that this EOS is quite different from the 
EOS for the pure hadronic phase as illustrated in Fig.3 in the case of an ideal pion gas. 

Thus, the lattice QCD data do really predict a minimum in the p/e- representation of 
the EOS whose position, according to Hung and Shuryak [0, defines the beam energy 
at which the fireball formed has the longest lifetime. This representation is significantly 
more sensitive, both to details of the lattice calculations and approximations involved 
therein than the conventional thermodynamic quantities s{T) and p{T). The simulation 
of the EOS with a crossover [|[] results in the softest point at £sp ~ 1.5 GeV/ fm^, which is 
noticeably higher than the lattice value Egp ~ 0.5 GeV/ fm?. Describing correctly the order 
and temperature of the deconfinement QCD transition for ub = 0, the mixed phase model 
1^, predicts Esp ~ 0.3 GeV/ frn? and the agreement with the lattice data may be further 
improved by a more accurate treatment of cluster-cluster interactions. All this implies 
that the proposed beam-energy tuning for identification of the deconfinement transition 

^These calculations within a hydrodynamic model are now in progress. 
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should be done at bombarding energies of Eiab ~ 3 — 5 GeV/A, much below the value 
Eiab — 30 GeV/A advocated by Hung and Shuryak [|l|. In the recent paper |jl6|, heavy 
ion collisions were considered within relativistic dissipative hydrodynamics with EOS 
taken from the lattice QCD calculations for baryonless matter. A strong enhancement 
of the width of the shock front has been found at Eiab — 6 GeV/A as a manifestation 
of the softest point. This value of the beam energy is in a fine agreement with our 



estimate, if one remembers that the old QCD data used in ref.[16| correspond to higher 
deconfinement temperature T^ec = 200 MeV. Note that these hydrodynamic calculations 
are complementary to the Hung and Shuryak estimate and show that the softest point 
effect influences not only the disassembly stage but also the compression stage. 

The mixed phase model predicts also a strong dependence of the EOS on the baryon 
density of the system: a minimum of the p/e function is washed out for ub ^ tiq. Since 
the state with e = Egp is a transitional one, we expect that the change of the fireball 
lifetime with Eiab will not be as large as predicted in Nevertheless, we would like 
to draw attention to heavy ion collisions at moderate {Eiab = 2 — 10 GeV/A) energies 
for studying the mixed phase of quarks and hadrons. As has been shown above, a pure 
hadronic EOS is quite different from the EOS predicted by the mixed phase model near 
the softest point. A possibility of forming the mixed quark-hadron state at energies 



Eiab ~ 2 — 10 GeV/A has been noted previously |T^. However, being based on the two- 
phase models, this consideration predicts a sharp decrease of temperature just above the 
deconfinement threshold resulting in the formation of a 'cold' plasma. Generally speaking, 
the known signals of quark-gluon plasma are not applicable to the cold plasma case. 
In contrast, the mixed phase model has no threshold and no temperature fall-off, thus 
such signals should persist but their strength will be proportional to an unbound quark 
abundance. It is worth mentioning that some enhancement of the A-hyperon production as 
a specific plasma formation signature has been observed at the energy as low as 3.5 GeV/A 



17|. Interferometry measurements deserve a special attention since they are sensitive to 



'granularity' of an emitting source |]rB|. In this respect, the correlation length or screening 
length of unbound quarks provides a new length scale additional to the source radius. 

We are indebted to W. Bauer, J. Cleymans, M. Gorenstein, E. Okonov, S. Pratt, 
F. Webber and G. Zinoviev for stimulating discussions and comments. V.D.T. acknowl- 
edges the warm hospitality of the theory group of GANIL, Caen where a part of this work 
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Figure captions 



Fig.l. Temperature dependence of energy density e and pressure p (relative to cor- 
responding values in the Stefan-Boltzmann limit) for SU (3) gauge theory of baryon-free 
matter with massive dynamical quarks. Curves are calculated within the mixed phase 



model P]. Triangles and squares are lattice data for Wilson and Kogut-Susskind |]l3l 
schemes of accounting for dynamical quarks, respectively. 

Fig. 2. The ratio of pressure and energy density p/e versus e. Notation is the same as 
in Fig.l. The dashed curve corresponds to the approximation used in ref.[||]. Data are 
given for e > 0.005 GeV fm~^ . 

Fig. 3. p/e-representation of the EOS for the baryonic matter predicted by the mixed 
phase model . Numbers near the curves show the baryon density in units of the normal 
nuclear density. The case of ideal pion gas is given by the dashed line. 
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